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| Abstract 

In this letter we establish a connection of Picard-type elliptic solutions of Painleve VI equation 
with the special solutions of the non-stationary Lame equation. The latter appeared in the study 
of the ground state properties of Baxter's solvable eight- vertex lattice model at a particular point, 
77 — 7r/3, of the disordered regime. 



^ ! 1 Introduction 

' The Painleve transcendents have numerous remarkable applications in the theory of integrable models 

£Sj ■ of statistical mechanics and quantum field theory (see, for instance, [1-3]). Mention, in particular, the 

calculation of the "supersymmetric index" and related problems of dilute polymers on a cylinder which 
lead to Painleve III [4,5]. These problems are connected with the finite volume massive sine-Gordon model 
with N = 2 supersymmetry. The lattice analog of this continuous quantum field theory corresponds to 
a special case of Baxter' famous solvable eight-vertex lattice model [6]. In this paper we continue our 
study [7] of this special model on a finite lattice and unravel its deep connections with Painleve VI theory. 

We consider the eight- vertex model on a square lattice with an odd number, N = 2n + 1, of columns 
and periodic boundary conditions. The eigenvalues of the row-to-row transfer matrix of the model, T(u), 
satisfy the TQ-equation [6], 



o 



- 1—1 

X 



T{u) Q{u) = <p(u - 77) Q(u + 277) + cf>(u + 77) Q(u - 2r)), (1) 



H 

d ' where u is the spectral parameter, 

cf>{u) =<(u|q), q = e i7TT , Imr > 0, (2) 

and #i(u|q) is the standard theta-function with the periods it and ttt (we follow the notation of [8]). 
Here we consider a special case 77 = tt/3, where the ground state eigenvalue has a very simple form [9, 10] 
for all (odd) N 

T{u) = <P{u), 77= J. (3) 

The equation with this eigenvalue, T(u), has two different solutions [11,12], Q±{u) = Q±(u,q,n), 
which are entire functions of the variable u and obey the following periodicity conditions [6, 13] 1 

Q±(« + 7r)=±(-l) B Q±(tt) J Q ± {u + ht) =q- N ' 2 e- iNu Q T (u), Q±(-«) = Q±(«). (4) 
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The above requirements uniquely determine Q±{u) to within a common it-independent normalization 
factor. It is convenient to rewrite the functional equation Q for Q±(u) in the form 



2lT 47T 

*±(u) + *±(« + y) + *±(« + y ! 



0, 



where 



*±(«) = *±(«,q,n) 



^" +1 (u|q) ( 



:(w,q,n), 



(•5) 



(6) 



z^(3u|q 3 ) 

are the meromorphic functions of the variable u for any fixed values of q and n. As shown in [7] these 
functions satisfy the non-stationary Lame equation 



6q — ¥(u,q,n) = { - ^ + 9n(n+l)P(3u|q 3 )+c(q,7i)}*(u,q,7i), 



(7) 



where the elliptic Weierstrass p- function, *P(v | e I7re ), has the periods it and 7re [8]. The constant c(q,n) 
appearing in J7J is totally controlled by the (u-independent) normalization of Q±(u). Note that this 
equation (in fact, a more general equation, usually called the non-autonomous Lame equation) arose 
previously [14-17] in different contexts. We will not explore these connections here. 
Let us quote other relevant results of [7]. Define the combinations 



Qi(«) = (Q+(u) + Q-(u))/2, Q 2 (u) = (Q+(u) - Q-{u))/2, 



such that 



l (u + 7r) = (-1)' 



(«)• 



(8) 



(9) 



Bearing in mind this simple relation we will only consider Qi(u), writing it as Q^(it) to indicate the 
n-dependence. Obviously the linear relations JSJ and Q remain unaffected if Q±(u) in 10 is replaced by 
Q^\u). The partial differential equation has, of course, many solutions. Here we are only interested 
in the very special solutions, relevant to our original problem of the eight-vertex model (by definition the 
functions Q^' \u) are entire quasi-periodic function of u with the periods 7r and 2txt). Introduce new 
variables 7 and x, instead of q and u, 



7 = 7(q) 



tfiU/3|q 1/2 ) 



x = 7 



tf 3 W2|q 1/2 ) 



, 2 



_i? 2 (V3|q 1/2 ) 
and new functions V n {x, z) instead of Q^(u), 

Qi n) (u) = W(q, n) ^{u/2 | q 1 / 2 ) ^ n {u/2 | q 1 / 2 ) V n (x, z 



^W2|q 1/2 ) 



(10) 



(11) 



where A/"(q,n) is an arbitrary normalization factor (which remains at our disposal). The properties of 
the functions V n (x,z) corresponding to the required solutions of (JSJ) and (0) can be summarized by the 
following 

Conjecture 1 

(a) The junctions V n (x, z) are polynomials in x, z of the degree n in x, 



V n (x,z) = J2r£\z)x k . 



fe=0 



while r^ n \z), i = 0, . .. ,n, are polynomials in z of the degree 

deg[ri n) (z)] < [n(n- l)/4 + fe/2j 

with positive integer coefficients. The normalization of V n {x, z) is fixed by the requirement 
Tn (0) = 1 and \x\ denotes the largest integer not exceeding x. 



(12) 



(13) 
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(b) The coefficients of the lowest and highest powers in x, corresponding to k — and k — n in 
read 

ri n \z) = r n (z, -1/3), r^(z) = r n+1 (z, 1/6) , (14) 

where the functions T n (z,£) (for each fixed value of the their second argument are determined by 
the recurrence relation 

2z(z - l)(9z - l) 2 [logr„(z)];' + 2(3z - l) 2 (9z - l)[log r„(z)]' z + 

2 T n+ l(z)T„_l(z) 



+8 



2n - 4£ - - 
^ 3 



-[12(3n - 6f - l)(n - 2$) + (9z - l)(n - l)(5n - 120] = 0, (15) 
luj'i/l i/ie initial condition 

t (z, = 1, n(2,0 = -44 + 5/3. (16) 
TTie functions T n (z, are polynomials in z for all n = 0, 1, 2, . . . , oo. 

The conjecture has been verified by an explicit calculation [7] of V n {x, z) for all n < 50. 

The polynomials 7-Vj(a;, z) can be effectively calculated using the algebraic form (see [7]) of the equation 
0. The first few of them read 2 

Vq(x,z) = 1, Vi(x,z)=x + 3, V 2 (x,z) = x 2 (l + z) + 5.x(l + 3z) + 10, (17) 

P 3 (x, z) = x 3 {\ + 3z + 4z 2 ) + 7ir 2 (l + 5z + 18z 2 ) + 7x(3 + 19z + 18z 2 ) + 35 + 21z. (18) 

As explained in [7], Eq.JJJ) leads to a descending recurrence relations for the coefficients in (|12J) . in the 
sense that each coefficient r^\z) with k < n can be recursively calculated in terms of rm\z), with 
m = k+ 1, . . . , n and, therefore, can be eventually expressed through the coefficient r^ (z) of the leading 
power of x. The conditions that this procedure truncates (and thus defines a polynomial, but not an 
infinite series in negative powers of x) completely determine the starting leading coefficient as a function 
of z. The above conjecture implies that these truncation conditions are equivalent to the particular case 
(£ = 1/6) of recurrence relation l)15[l. (|16|l . A similar statement applies to the (smallest power in x) 
coefficient r^ of V n (x, z). The results quoted above are due to [7] (except for the expression for r^ in 
(J14JI , which is new). 

In this paper we show that the recurrence relation H15[) exactly coincides with that for the tau- 
functions associated with special elliptic solutions of the Painleve VI equation, revealing hitherto unknown 
connections of the eight-vertex model and the non-stationary Lame equation Q to Painleve VI theory. 

Before concluding this introduction, let us briefly mention some problems related to the special eight- 
vertex model considered here. In the trigonometric limit (q — > 0) this model reduces to the special 
six- vertex model (with the parameter A = —1/2), which is closely connected with various interesting 
combinatorial problems [18-21], particularly, with the theory of alternating sign matrices. The fact 
that the polynomials V n (x,z) have positive integer coefficients makes it plausible to suggest that these 
coefficients have some (yet unknown) combinatorial interpretation. 

In the scaling limit 

n -> oo, q -> 0, t = 8 riq 3/2 = fixed, (19) 

the functions © essentially reduce to the ground state eigenvalues Q±(0) = Q±(0,t) of the Q-operators 
[11,22] of the restricted massive sine-Gordon model (at the so-called, super-symmetric point) on a cylinder 
of the spatial circumference R, where t = MR and M is the soliton mass and the variable 9 is defined 
as u — ttt/2 — i0/3. With a suitable ^-dependent normalization of Q±(9), Eq.0 then reduces to the 
"non-stationary Mathieu equation" [7] 

t§- t Q±(0,i) = {^-^ 2 (cosh20-l)} e±(0,t). (20) 



2 The higher polynomials with n < 50 are available in electronic form upon request to the authors. 
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This equation determines the asymptotic behavior of Q±{0) at large 9 



log Q±{6) = -jte 6 + \ogV ± (t) + 2 (d t \ogV ± {t) - t/8) e' + 0{e- w ), 6 -» +oo, (21) 

where D±(i) are the Fredholm determinants, which previously appeared in connection with the "su- 
persymmetric index" and the dilute polymers on a cylinder [4,5,23,24]. Note, in particular, that the 
quantity 

F(t) = ±U®, ^)=log^|, (22) 
describes the free energy of a single incontractible polymer loop and satisfies the Painleve III equation [23] 



1 d d , s 1 

-— t— U(t) = - smh2U(t). 

tdt dt v ' 2 w 

2 Painleve VI equation 

The Painleve Py/(a, 0, 7, S) is the following second order differential equation [25,26] 



l"(t) 



1 / 1 



2\q(t) q(t)-l q(t)-t 
q(t)(q(t)-l)(q(t)-t) 



t 2 {t-l) 



a + (3 



t 



t (t-1) q(t) - t 



<?'(*)- 



t- 1 



q(t) 2 



+ 5- 



t(t-i) 



{q{t)-lY (q(t)-t) 



where the following parameterizations of four constants [27] is chosen 



1 







1 



1 



,2 



1 



(1 



2 00 ' ' 2 u ' ' 2 ^ 2 
Ko = bi+b 2 , n\ = bi - b 2 , K oc =b 3 -b i , 9 = 63 + 64 + 1. 
This equation is equivalent to the Hamiltonian system Hvi(t; q,p) described by the equations 

dH 
dq 

with the Hamiltonian function 

1 



dq dH dp 
dt dp ' dt 



Hvi{t]q,p) 



q{q-l){q-t)p A 



t(t - 1) 

-{k (q -l)(q-t) + /si g(g - t) + (0 - l)q(q - l)}p + K (q - t)}, 



where q = q(t) , p = p{t) and 



k = i(/t + ki + - 1) 2 - i«4> 



One can introduce an auxiliary Hamiltonian h(t), 



h(t) = t(t - l)H{t) + e 2 (6i, 63, 64) * - ^ 2 (b u b 2: b 3 , b^ 



(23) 



(24) 

(25) 
(26) 

(27) 



(28) 
(29) 

(30) 



where ei(xx, . . . ,x n ) is the i-th elementary symmetric function in n variables. 

Okamoto [27] showed that for each pair {q(t),p(t)} satisfying l|27|) . the function h(t) solves the ~Eyi 
equation, 



h'(t) t(l - t)h"(t) + h'{t)[2h(t) - (2t - l)h'(t)} + 
and q(t) solves Py/(a, (3, 7, 5), J2U- 



n(v(i) 



(31) 
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Conversely, for each solution h(t) of l|31l) . such that 4sh(t) ^ 0, there exists a solution {q(t),p(t)} of 
(|27() . where g(i) solves (|24ll . An explicit correspondence between three sets {<?(£), q'(t)}, {q(t),p(t)} and 
{/i(i), /i'(t), h"(t)} is given by birational transformations, which can be found in [27]. 

The equation ~Evi is a particular case of a more general class of the second order second degree 
equation 

y"(t) 2 = F(t,y(t),y'(t)) (32) 

where F is rational in y(t), y'{t), locally analytic in t, with the property that the only movable singularities 
of y(t) are poles. Such equations were classified in [28], where Ea. (l31|l is refered to the "SD-I type". 
Originally, we have obtained equations of the form Il-'i2[l for asymptotics of the polynomals "P n (x, z) and 
only then reduced them to (|31|1 and Ij24(l . 

The group of Backlund transformations of Pyi is isomorphic to the affine Weyl group of the type D4. 
It contains the following transformations of parameters (only 5 of them are independent) 

wi : b\ <-> &2, w 2 ■ b 2 <-> b 3 , w 3 : b 3 <-> 64, w 4 : b 3 — > —63, 64 — > -fe 4 , (33) 

x 1 : Ko <-* Ki, x 2 : k <-> k^, x 3 : n a ^ 9 (34) 
and the parallel transformation 

Z 3 : b = (h, b 2 , 63, 64) -> b+ = (61, 63, 63 + 1, 64). (35) 

Here we shall use only two transformations: x 2 and Z3. The canonical transformation x 1 corresponds to 
a simple change of variables [27] , 

xl : (q,p,H,t) -» (q-\eq-q 2 p,-H/t 2 ,l/t) 7 e = ]^{n + n x + 9 - I + k^) . (36) 
The birational canonical transformation, 

{<?,p} = {?(b),p(b)} - {<Z+,p+} = { 9 (b+),p(b+)} , (37) 
corresponding to ^3, can be obtained from the observation that [27] 

h+(t, q + ,p+) = h{t, q,p) - q{q - l)p + (6j + b A )q - -(b x + b 2 + b 4 ). (38) 
Here we shall only give the expression for q in terms of {q + ,p + } 

q{p 2 + q+(q+ - l)(q+ - 1) + +p+[2(l + 61 + b 3 )q+(l - q+) + (61 + b 2 )(q+ - t) + 2b 1 (t - l)q+] + 

+ (1 + h + 63) [(1 + 61 + b 3 )(q+ - 1) + t(l + b 3 - b x ) + 61 - 62] }- 

-<[p+(g+ - 1) - 1 - 61 - 63] [p+<z+(<z+ - 1) - (1 + h + b 3 )q+ + b! + 62] = 0. (39) 
Later we will need the transformation x 2 l 3 x 2 

x 2 l 3 x 2 : {b} = {61, 6 2 , 63, 64} -> {b} = {6j + \M + \ b 3 + \M + \}- (40) 
Combining (|36I39|I one obtains 

Q = {p 2 (v - *)(? - 1) +P ( 2 + h + h + b 3 + bi + (h - b 2 )t - (2 + 26 x + 6 3 + b 4 )q) 
+ (l + b 1 +b 3 )(l + b 1 +b i )}/{(l + b 1 +b 3 +p{l-q)){l + b 1 +b i +p{l-q))} , (41) 

where q = q(h) and q = q(b), p = p(h). 
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3 Special elliptic solutions 

It goes back to Picard [29] that if parameters l|26|) satisfy 

h = b 2 = 0, 6 3 = b i = -1/2, (42) 
then a general solution of Pvi reads (see also [30]) 

q{t) = p{c!bJi + C2LJ 2 ;U)1,0J2) + (43) 

where p(it; wi,o>2) is the Weierstrass elliptic function with half-periods 1^1,2, cj.,2 are arbitrary constants 
and a;i i2 (i) are two linearly independent solutions of the hypergeometric equations 

t(l - t)u"(t) + (1 - 2t)u'{t) ~ -u(t) = 0. (44) 

It is convenient to choose 

wi(t) = ~ 2F 1 (\, \; l;t) = K(^ ); ^ {t) = f | 1 , i; 1; 1 - t) = iK'^/S), (45) 

where K and K' are the complete elliptic integrals of the modulus k = t x l 2 . 
Using expressions for invariants of the Weierstrass function e\, e%, e% 

ei = l~^p, e 2 =f-^yl, e 3 = (46) 
we can rewrite Picard's solution of Pvi as 

g(t) = ns 2 (ciK + c 2 iK',A:), k = t 1/2 . (47) 

Let us choose 

ci = 1, ca = - (48) 

and denote corresponding solution of (|24|1 as qi(t). Using addition theorems for elliptic functions it is 
easy to show that qi (t) satisfies the algebraic equation 

q$(t) -*tq\{t) + Uq\{t) -4iq x (t) + t 2 = 0. (49) 

In fact, this algebraic solution of Pvi is very special. It solves <|24[) not only for the Picard's choice of 
parameters (|4*2*)l but also for 

&i=£-J> ^2 = 0, h = i-\, b i = 2£-^ (50) 
6 3 6 

where £ is an arbitrary parameter. This happens because Pvi, l|24|l . splits up into two different equations 
which are both satisfied by I^J. 

Using above formulas we can easily find expressions for p±(t) and h±(t) corresponding (|49|l 

(t) _ (l-3Q(3t-2t qi (t)-g 2 (t)) l-2t (l-2Q(l-3Q(t + g 1 (t)-2t gi (t)) 

Pl[ ' 6t{ qi (t)-l) 2 lU 72 4 ' 

Now we shall assume that this solution corresponds to the case n = 1 and apply Backlund transformation 
x 2 (lz) l ~ n x 2 to obtain a series of solutions {q n (t),p n {t),hn(t)} with parameters 

1 n , 1 n 1 n , l^n , 

&1 =3 + ^2< 62 =2-2< 63 = -6 + ^2< 5 ^-3 +2 ^2- < 52 > 

At this stage we are ready to establish a connection with the elliptic parametrization from the first section. 
Let us assume that the elliptic nome q in (J2J) , Ij6l7(l and llOjl is given by 

f 2K'(fc)i ,/n 
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where K(fc) and K'(fc) are denned by l|45|l . 

Using Landen transformation for elliptic functions it is easy to obtain the following rational parametriza- 
tion for z = l/j(q) 2 defined in JTUJ), fTTft. and for t, qi(t) in (gSJ 



'7 



(l-7)(3 + 7) 3 
(l + 7)(3-7) 3 ' 



Qi(t) 



(l-7)(3 + 7) 
(l + 7)(3-7) : 



(54) 



7 2 (q) (3-7)7' 
in terms of a new parameter 

7^7(q 1/2 ), (55) 

defined by l|10[l with q replaced by q 1 / 2 . Note that such parameterization of Picard's solutions of Pyi 
with the above choice (j48(l of the parameters ci and C2 has already appeared in [30,31]. 
From these formulas one can get an explicit connection of variables t and z 



t = 



(z-l)(l-9z) 3 
32~2 



27 'z 1 - 182 - 1 



v/(l-z)(l-9z) 3 



(56) 



Now we can construct a sequence of r-functions associated with a series of auxiliary Hamiltonians 
h n (t). It appears that corresponding r-functions are polynomials in variable z. 

First, let us introduce a sequence of functions cr n (z) considering them as functions of z 



an{z) = r z y — 



h n (t) + ^(2t - 1) + (n - l) 2 



t-1 l-9z 



(57) 



+("-!)(*- 12) 



1-t- 



t(l - 3z) 



1 - * 



v/(l-z)(l-9z) 
Comparing it with (|51|l and using l|54|l it is easy to see that 

(7i(z) = 0. 

Then using Backlund transformations 2 and x 2 l^x 2 it is not difficult to show that 

at(z) = 0, i = 0,l,2. 



(58) 



(59) 



The easiest way to do that is to calculate hi(t) , i = 0,2 in terms of h\(t) (|5ip. substitute into (|57|l and 
use a rational parametrization Q54)). 

Now let us introduce a family of r-functions r„(z,^) via 



<J n (z) = — [log r„ 



and fix a normalization for n = 0, 1, 2 as 



-4^ + 5/3, r 2 (*,0 = 3(2£-l)(3£-l). 



(60) 



(61) 



Using Okamoto's Toda-recursion relations for r-functions for Pyi, generated via successive applications 
of parallel transformation Z3 [27], one can show the recurrence relation for r„(z, £) exactly coincides with 
(|15(l . Thus, we showed that the leading coefficient and the constant term of P n (x,z) (considered as 
polynomials in x) can be expressed in terms r-functions for special solutions of Pyj. 

At the moment we do not have a complete proof of the polynomiality of r„(z). Note, that this property 
takes place provided that two successive r-functions r„(z) and r n +i(z) do not have a nontrivial common 
divisor (which is a polynomial in z). 

One of the challenging problems is to find a determinant representation for r„(z,£) similar to that 
known other polynomial solutions of the Painleve equations. It could help to clarify the structure of 
V n (x, z) and possibly to establish a connection with problems of combinatorics. 
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